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We investigate HBT radii based on the numerical solutions of the hydrodynamical 
model which are so tuned as to reproduce the recent experimental data at the CERN 
SPS. Comparing the sizes of freeze-out hypersurface with HBT radii, we discuss dy- 
namical effect of the systematic flow on the apparent HBT radii. Finally we compare 
HBT radii of the QGP phase transition model with those of the hadron gas model 
without phase transition. 

Hanbury-Brown Twiss (HBT) effect [0] is a well known quantum effect which enables us 
to estimate the source size through the two-body correlation of the emitted particles. In the 
high energy nuclear collision, the correlation experiments of pions, kaons, etc are promising 
for obtaining the knowledge of the space-time size of the hot and dense region during the 
reactions. However, in the nuclear collisions, the reaction takes place highly dynamically and 
the particle source is not static. The produced hot fire-ball is expected to expand rapidly 
and to cool down in very short time period. Hence, the meaning of the "size" of the source 
(fire-ball) is not trivial. 

In previous papers, we analyzed numerically (3+l)-dimensional hydrodynamical model 
f§] and applied it for the recent experimental data of of CERN WA80. In this paper, based on 
the numerical solutions we evaluate two-particle correlation from which we estimate the HBT 
radii. Then we compare these results with the size of freeze-out hypersurface. In particular, 
we analyze the effect of systematic flow of the fluid. Finally, in order to investigate the 
signature of QGP, we compare the HBT radii of the phase transition model and those of the 
hadron gas model without phase transition. 

An annihilation operator of the particle (e.g. pion) emitted from a source J(x) is given 

as 
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where k is momentum of the particle and is on-shell frequency, = y k 2 + m 2 . The 
subscripts 'out' and 'in' correspond to out-going field and in-coming field, respectively. Here 
we assume J(x) is a c-number source with random phase, i.e., 

J(i)=7ge^ ) , (2) 

where (p(x) is a random number. The statistical properties of 4>{x) are given by 

< e i*(*) e -i<M >= s\x-y), (3) 

and Gaussian reduction holds for the higher order correlations. Then one-particle spectrum 
is given as, 

W{1 ^ = ^k = < { < -l a OUt( fc ) a OUt( fc )|0m >} > 



where < ... > is the average over the random phase and |0j n > is the vacuum state of the in- 
coming field. The c-number source J(x) is defined from consistency in the hydrodynamical 
model, i.e., one-particle spectrum should be given by the thermal distribution, 

dN r U^k 



d 3 k JT=Tf /i (27r) 3 a; fc ' 

where is the four-velocity of the fluid and f(E,T) is Bose-Einstein distribution function. 
Comparing the above formula with Eq. 01, we can assign c-number source function as, 



J(x) = pU^fiU^T) 



(6) 

T=T t 



The two-particle distribution is given by, 

W(k u k 2 ) =< {< m |aJ )ut (fc 1 )aJ )ut (fc 2 )a out (fc 1 )a out (fc 2 )|0 i „ >} >, (7) 

and the correlation function which is comparable to experimental data is obtained after 
integrating about the average momenta, K = kl + k2 ; as 

w jd?Kw{k 1 )w{k 2 y y ' 

Following the conventional manner, we define apparent HBT radii, R s [d e , R ut an d -^long 
through the Gaussian fitting 0], 

Cfit(q) = 1 + Aexp j-- (i?g ide g s 2 ide + ^out^out + ^longSfong) } > ( 9 ) 




Figure 1: Example of the source distributions as a function of t, z and r, respectively. The 
solid line stands for the phase transition model and the dashed line stands for the hadron 
gas model. 



where A is chaoticity and q s i^ e , (fout an d Qlong are the side- ward component, the out- 
ward component and the longitudinal component of q, respectively. Assuming cylindrical 
symmetry of the Gaussian source function, and taking the limit of K ^> q, we have HBT 
radii as follows, 

^side = Ar 2 , ( 10 ) 
R 2 ont = Ar 2 + f3 2 At 2 , (11) 
4>ng = ^ 2 + p 2 L At 2 , (12) 

where Ar, Az and At are the transverse width, the longitudinal width and the temporal 

dE TC 

width of the source, respectively. (3t and (3l are the velocities as usual, (3t = and 

Pl — qx ■ Though the freeze-out hypersurfaces of our numerical solutions are not the 
Gaussian shape (Fig. |I|), we evaluate Ar, Az and At as the deviations of the volume element 
on the freeze-out hypersurface. 

Figures [I] and [| show HBT radii as a function of the hypersurface size. In order to 
clarify the effect of the hydrodynamical flow, in these figures we compared the two sets of 
data: one is the HBT radii as functions of the hypersurface size and the other set is the 
results of static picture. Data of static picture are evaluated with putting all four velocity as 
f/ M = (1,0,0,0) by hand. Four plots in these figures correspond to the S+S, S+Cu, S+Ag 
and S+Au collisions, from left to right, respectively. 

In figure ||, R s \& e is proportional to Ar in both data, hence, we may say that R s \<\ e can 
indicate Ar correctly even if the systematic flow exist. Because of the non-Gaussian shape 
of the source, the value of R s \<\ e doesn't equal to the Ar. On the other hand, though i? Q ut 
is proportional to J Ar 2 + fi\ At 2 (Fig. without systematic flow, R ou t, in the case with 
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Figure 2: R S ide as a function of Ar. 
Squares stand for the source without 
flow and crosses stand for the one 
with flow. 



Figure 3: 

Rout 

function of 

(Ar 2 + /3j.At 2 )a. Squares stand for 
the source without flow and crosses 
stand for the one with flow. 



flow, is not proportional to yAr 2 + /3^At 2 and becomes smaller than it. The systematic 
flow causes the reduction of apparent out-word HBT radius, i? Q ut- When flow exists, -R} on g 
also become smaller. However R\ on g is not proportional to the size of source even if flow 
doesn't exist. 

We also calculate the two-particle correlation numerically with use of the hydrodynamical 
model of the hadron gas model without phase transition J3J. Figures |^, ^| and [7] show 
the comparison between the QGP-phase transition model and hadron gas model. In these 
figures, four plots again correspond to the S+S, S+Cu, S+Ag and S+Au collisions, from left 
to right, respectively. The hadron gas model gives larger -Rgj^g than the phase transition 
model, because Ar of the hadron gas model is larger than that of the phase transition model 
as shown in Fig. [I]. In the cases of R on t and R\ on g, the results of these models seem to be 
hardly distinguishable. 

In this paper, we compared HBT radii with the size of source calculated numerically, 
-fi'side corresponds to the size of source. However, i? ou t and R\ on g become smaller than the 
size of source because of the systematic flow. We also compared HBT radii of the phase 
transition model with those of the hadron gas model. According to our numerical results, it 
seems very difficult to distinguish these models with HBT radii only. 

The authors are indebted to Prof. I. Ohba and Prof. H. Nakazato for their fruitful com- 
ments. They also thank to many discussions with members of Waseda Univ. High Energy 
Physics Group. 
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Figure 4: Riong as a function of 
(Az 2 + PlAt 2 )^. Squares stand for 
the source without flow and crosses 
stand for the one with flow. 



7.5 
7 

6.5 
6 

5.5 
5 



20 



+ 
□ 



+ 
□ 



110 

Target Mass Number A 



200 



Figure 5: R S ide as a function of target 
mass number. Crosses stand for the 
phase transition model and squares 
stand for the hadron gas model. 
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Figure 6: R out as a function of target 
mass number. Crosses stand for the 
phase transition model and squares 
stand for the hadron gas model. 



Figure 7: Ri ong as a function of target 
mass number. Crosses stand for the 
phase transition model and squares 
stand for the hadron gas model. 
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